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ABSTRACT 

oo 

^ I We examine the effect of mass and force resolution on a specific star formation 

i_^ _ (SF) recipe using a set of A'^-body/Smooth Particle Hydrodynamic simulations of 

O I isolated galaxies. Our simulations span halo masses from 10^ to 10^^ Mq, more 

' than four orders of magnitude in mass resolution, and two orders of magnitude 

in the gravitational softening length, e, representing the force resolution. We 
is^ ^ examine the total global star formation rate, the star formation history, and 

^ I the quantity of stellar feedback and compare the disk structure of the galaxies. 

Based on our analysis, we recommend using at least 10^ particles each for the dark 
matter and gas component and a force resolution of e ^ lO'^R^r when studying 
global SF and feedback. When the spatial distribution of stars is important, the 
number of gas and dark matter particles must be increased to at least 10^ of each. 
Low mass resolution simulations with fixed softening lengths show particularly 
weak stellar disks due to two-body heating. While decreasing spatial resolution 
in low mass resolution simulations hmits two-body effects, density and potential 
gradients cannot be sustained. Regardless of the softening, low-mass resolution 
simulations contain fewer high density regions where SF may occur. Galaxies 
of approximately 10^° Mq display unique sensitivity to both mass and force 



c^ 
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resolution. This mass of galaxy has a shallow potential and is on the verge of 
forming a disk. The combination of these factors give this galaxy the potential 
for strong gas outflows driven by supernova feedback and make it particularly 
sensitive to any changes to the simulation parameters. 

Subject headings: hydrodynamics, stars:formation, methods: numerical, galaxies: 
structure, galaxies: evolution 



1. Introduction 



The formation of stars is a highly non-linear process that is dependent on a number of 
different physical m echanisms, po ssibly including but riot lirnited to: gravitational (IZinnecker 
1984 : iLarsonlllQSSl ) or turbule nt (iPadoan fc Nordlundlbooa ISate fc Bonnellll2004h fr agmen- 
tation, co mpetitive accretion (Bonnell fc Batell2006l) . magnetic fields (IShu et al.|l2004), stellar 
feedb ack (ISilklll995l : lAdams fc Fatuzzdl 19961 ). and the coalescence of protostars (IMurrav fc Lin 
1996 ) . On larger sca les, SF is driven by gravit ational interactions betw een galaxies (JBarton et al 



2000; Liet a 



2008 ), external gas accretion ( iKauffmann et alj|2006[). and disk instabilities 



such as bars (jHo et al.lll997l : ISheth et al.ll2005[ ) and spiral arms ( lRobertslll969l : lTamburro et al. 



20081). Furthermore , energy from active galactic nuclei and supernova feedback (JMartiru 



1999 



Kaviraj et al.l 120071 : iBundy et al.l 120081 ) regulate SF. Hence, SF is determined by the inter- 
action between many large and small scale physical processes. 

SF is commonly studied on a galactic scale using A'^-body/Smoothed Particle Hydro- 
dynamic (SPH) simulations because of their ability to model processes over orders of mag- 
nitude in space and mass while automatically resolving denser areas to a higher degree. In 
fact, it is only very recently that an Eulerian hydrodynamic simulati on has followed th e 



formation of a galaxy in a cosr nological context to a redshift of zero (Gibson et al. 



Brooks et al 



2009). 



20091 ) 



Cosmological SPH simulations ( iKeres et al.ll2005l : iGovernato et al.ll2007l : 
contain one or more galaxies forming within a larger, cosmologically determined dark mat- 
ter (DM) structure an d are crucial for replicating the tidal interac t ions with other galaxie s 
that may trigger SF ( lOkamoto et al.l l2005t iKaufmann et al.l l2006t iGovernato et al.ll2007l ). 
These cosmological sir nulations are used for studying the connection between SF and UV 
background radiati on ( jBuUock et al.ll2000l : lKravtsov et al.ll2004j: JMoore et al.ll2006l) . mergers 
( JAbadi et al.ll2003), tidal and ram pressure strip ping ( JMayer et al.ll2006l ). and gas accretion 
by galaxies ( jKeres et al.ll2005t JBrooks et al.ll2009l ). On a smaller scale, simulations of isolated 
galax ies give a resolved naodel of SF within a galaxy and its effect on the galactic struc- 



2007 



Roskar et al.ll2008[ ). 



ture fIValenzuela &: Klvpinll2003l : JDebattista et al.ll2004J : JTasker fc Brvanll2006 



Stinson et al. 
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Even in high resolution SPH simulations of the smallest isolated galaxies, the formation 
of individual stars is unresolved. In order to describe the non-linear, sub-resolution process 
in SPH simulations, SF is implemented using sub-grid models wherein the complex physical 
processes are averaged over the scale of the smallest resolution element. Comparisons with 
Eulerian grid models of SF in th e interstellar medium d emonstrate the viability of sub-grid 
SF methods in SPH simulations (JTasker fc Bryan! |2006| ). Given basic information about the 
surrounding gas, such as density and temperature, statistical models of SF can reproduce 
most observed global SF properties of galaxies in well-resolved simulations. The question 
remains, however, what is necessary for a simulation to be "well-resolved"? Although the 
qualitative behavior of SF in SPH simulations has been discussed in previous papers, here we 
will establish the minimum spatial and mass resolution necessary to m odel different aspects 
of SF in galaxies simulated with the widely used code GASOLINE ( IWadsley et al.l |2004[ ) 
over a wide range of masses. As well as delineating the behavior of SF with GASOLINE as 
a function of resolution, these results may be applicable to other SPH simulations in which 
individual gas particles represent single-phases of the ISM, SF occurs in suitably cold and 
dense gas particl es , and energy from sup e rnova is transferred to n e arby gas part i cles, e .g. 
Katz et all Jl996h : Ixhacker fc CouchmanI J200lh : IXeres et all ^W^ : ISaitoh et aP J2008[ l. 



Stinson et al.l ( 120061 ) . hereafter S06, described a specific sub-grid model for SF and a 



blastwave model for stellar feedback in GASOLINE, laid out a physical justification for the 
models, and presented several different tests of the SF and feedback models. To judge the 
effects of different SF criteria and model parameters, S06 used an isolated Milky Way-mass 
disk galaxy. When using the best-f it parameters, t his model galaxy was able to repro- 
duce the Schmidt law, T^sfr oc ^I'as ( iKennicuttI Il998l ) . with a near constant star formation 
rate (SFR). S06 further simulated this isolated galaxy at three different mass resolutions 
(2000, 10000 and 50000 initial gas particles and a gas disk mass of IO^Mq). While fewer 
stars formed in simulations run at lower reso l ution , the global SFR started to converge 
for 10000 initial gas particles. iBrooks et al.l (l2007l ) found that the total amount of SF 
and the metallicity converged for oyer a fe w thousand DM particles per galaxy in cosmo- 
logical simulations. iGovernato et al.l (|2007l ) further tested this SF recipe when calibrating 
the SF parameters through comparisons of the SFR and velocity dispersion of an isolated 
dwarf and a MW-mass galaxy to observations. This SF and feedback recipe, albeit using 
a different set of param eters, has been successful in reproduci ng the TuUy-Fisher Relation 
(IGovernato et al.ll2007n. the internal s tructure of disk galaxies (IGovernato et al.ll2004i l2007l : 



Roskar et al. 



12008 



Stinson et a 



a systems ([Pontzen et al. 



(IBrooks et al.ll2007n . 



2008 



20091) . the column density distributions of damped Layman 



2010l ). and the observed stellar mass-metallicity relationship 



Despite this success, it is possible that SF in galaxy simulations outside of the previously 
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tested range of resolution could behave inconsistently. ISaitoh et al.l (|2008[ ) found that when 
using a similar SPH/N-body code, in very high resolution simulations (3.5 x 10^ — IO^Mq 
per particle), increasing the threshold density for SF and decreasing the efficiency resulted 
in more realistic disks which were both thinner and clumpier. If simulations with different 
resolutions show different properties, such as a range of gas densities, the SF recipe may need 
to be adjusted at these extremes. As well as affecting the robustness of the particular sub-grid 
SF model used, resolution also affects the dynamics of galaxy s imulations and impacts the 



SF through the alteration of the density structure of the galaxy (ISpringel fc Hernguistl 12002 



Governato et al.ll2007l: iKaufmann et al. 



2007 



Naab et al.ll2007l : iFovle et al 



20081). Unlike the 



recipe put forth in ISpringel fc Hernguistl ( l2003a[ ) , which includes a sub-resolution multiphase 
model of the ISM, the S06 recipe assumes a singl e pressure, density and temperature for the 
entire gas particle. ISpringel fc Hernguistl (l2003al ) developed their multiphase model in part 
to minimize resolution effects but at the cost of small scale spatial structure of the multi- 
phase ISM. Because GASOLINE uses a single phase ISM for each particle it is especially 
important that the density structure of the ISM be well resolved to correctly model the gas 
associated with SF. Understanding the impact of resolution becomes particularly important 
when analyzing large-volume, cosmological simulations that contain a number of galaxies 
with different masses, each resolved to different degrees. The wide range of galactic masses 
in these cosmological simulations mean SF must remain accurate in a variety of different 
galactic environments, from large galaxies with high densities of gas to small ones from 
which gas is easily expelled. It is, therefore, necessary to establish how SF behaves at 
different mass and spatial resolutions over a range of galactic masses. 

There is a large body of work on the effect of both mass resolution (the mass of the par- 
ticles and inversely related to the number of particles) and force resolution (the gravitational 
softening length, e) on the dynamics of SPH simulated galaxies. The effects of insufficient res- 



olution can be largely group ed as those relating to over-cooling (jHutchings fc Thomas 



2007 



oiution can De largely group ea as tnose reiatmg to over-cooimg i|iiutcnmgs 
Springel fc Hernguistl 120021) . two-body heatin g (Steinmetz fc Whitelll997l: 



Maver et allboosh . artificial vi scositv kaufmann et aP I2OO7I: IWadslev et aP boosh 



2000 



Governato et al. 



and small-scale gravitational support ( JBate &: Burkertlll997l ). The first two effects primarily 
relate to the mass resolution whereas the third is connected to the force resolution. These 
effects can make the galaxy either more or less conducive to SF and they can even impact 
each other. Understanding how each alone impacts the dynamics of a simulated galaxy, 
however, is the first step to understanding how together they may infiuence SF. 

The first potential problem for accurately modeling SF in low resolution SPH simula- 
tions is that over-cooling of gas in galactic halos simulated with small numbers of particles 
( ISpringel &: Hernguistl 120021 ) may allow the gas to re ach the density threshold fo r SF too 
guickly This cooling, also known as in-shock cooling ( iHutchings fc Thomasll2000l ) happens 
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when low resolution causes the broadening of shock fronts and contact discontinuities and, 
as a result, the gas moving through them is heated relatively slowly. This longer heating 
time-scale gives the gas the opportunity to cool radiatively, which reduces the post-shock 
cooling time by as much as a factor o f two and lowers the maxira um shock temperature by 
up to 50% in shock-tube simulations (JHutchings fc Thomaj 120001 ) . 



In low mass resolution simulations, two-body heating and the processes related to it 
may affect SF by changing both the dynamics and the gas physics. Two-body heating 
results from approximating a smooth potential with discrete masses and causes an arti- 
ficial increase of kinetic energy. In particular, the presence of large collisionless bodies, 
such as dark matter (D M) particles, can heat the disk and destroy galactic disk structure 
(JLacey fc Ostrikerlll985[ ). Even when the rotational disk is not fully degraded through two- 
body heating, numerica l relaxation can transfer angular momentum f rom the disk to the halo 
(JKaufmann et al.l 120071 ). This result s in a weaker disk corn ponent ( iGovernato et al.ll2004j ) 
and excessively concentrated disks (IGovernato et al.l 120071 ). In addition to changing the 
dyna mics, two-body heating in creases the gas temperature through artificial viscous dissipa- 
tion (jSteinmetz fc Whitelll997l ). making it increasingly difficult for gas to reach temperatures 
below the temperature threshold for SF. 

One final way in which the small number of particles in a simulation may affect SF is 
through the reduced spatial resolution of the gas. Simulations with small numbers of particles 
have their gas properties calculat ed over larger snaoothi ng lengths, h, which in turn adds 
artificial viscosity and dissipation (JKaufmann et al.l 120071 ) and results in an increased loss of 
angular momentum and smaller, more concentrated disks. The smoothing of gas properties 
over larger volumes also results in fewer regions with sufficiently cool temperatures and high 
densities for SF. 

The ability of regions to collapse to form stars in A'^-body simulations is likewise sensitive 
to the force resolution, i.e. the gravitational softening length. Large values of e spread the 
potential ove r a larger volume, resu lting in decreased density fiuctuations and producing 
puffier disks (JKaufmann et al.l 120071 ). Increasing e, therefore, decreases the amount of gas 
able to reach sufficient densities for SF. Larger values of e also hinder the formation of bars 
as bars are only formed in the non-axisymmetric modes of the density field that are larger 
than e. Without bars, spiral stru cture is damped and less angular momentum is transferred 
outward (JKaufmann et al.ll2007l ) causing a lower concentration. Without the high densities 
of gas produced in bars, spiral arms, and in the center of the galaxy, SF may be reduced. 

When the Jeans length is comparable to either e or h the relationship between them 
becomes particularly important. When e ^ h, gravitational forces on scales smaller than 
the softening length are reduced, causing artificial support. Conversely, when e < h gravita- 
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tional forces dominate over pressure forces (JBate fc BurkertI 119971 ) . Reducing e relative to h 
increases the relative gravitational importance, while increasing it results in greater pressure 
support. Once the Jeans length begins to be unresolve d, gas dynamics, and through that 
SF, become dependent on the relative values of e and h ( iBate &: BurkertI 119971 ). The extent 
to which relative values of e and h affect SF is itself dependent on the implementation of the 
SF recipe. For this reason, it is important that the Jeans mass and length be greater than 
the respective mass and softening lengths of the star forming particles. 

As well as the effects of resolution on SF, the effects on the supernova feedback that 
regulates it must also be understood. All SPH feedback recipes must contend with the feed- 
back energy transferred to the nearby dense gas being radiated away too quickly, resulting 
in little hydrodynamical effect (JKatz et al.lll996[ ). One commonly used solution is to inject 
some of the SN energy in the form of kinetic rather than thermal energy to "wind particles" 
( ISpringel fc Hernquistll2003al ). Energy from feedback not used in the wind is necessary to 
maintain the equation of state in the presence of radiative losses. The effect of mass reso- 
l ution on the wind produced in this recipe have been examined by iDalla Vecchia fc Schaye 
( 120081 ) ■ among others. A second solution to simulating the sub-grid adiabatic phase of su- 
pernova explosions is to temporarily di sabling ga s cooli ng. This is the approach taken by 
the blastwave model, first presented in iGerritsenI (119971 ) and iThacker &: CouchmanI (J2000[ ) 
and later developed for Gasoline in S06. As no kinetic energy is explicitly transferred in this 
recipe, resolution will have a different effect. For example, in low mass resolution simulations 
using the S06 recipe, supernova events occur less frequently but are more powerful. In such 
simulations, the energy released per event is greater, that energy is distributed over fewer 
particles, and the cooling for each particle is disabled for a greater period of time (S06). 
While this recipe produces the same galaxy-and-time-averaged quantities irrespective of res- 
olution, (for instance, cooling is on average disabled for the same total mass of particles in 
a high resolution simulation as in a lower resolution one at any given time), it is an open 
question as to whether the different spatial and temporal distributions of supernova events 
in low and high mass resolution simulations will produce structurally different galaxies. Fur- 
thermore, the high spatial detail produced by using large numbers of particles may allow for 
the formation of non-spherically symmetric structures, such as galactic winds. Such struc- 
tures would l ikewise be affected by resolution induced changes to the artificial viscosity of 
gas particle s (IKaufmann et al.ll2007r) and changes to the disk scale height and the geometry 



of the disk ( Mac Low fc McCrav 



19881 : iMac Low et"aDll989h . 



In this paper, we explore the effects of force and mass resolution on S06 SF and stellar 
feedback in different mass galaxies. When analyzing galaxies with this SF recipe, it is 
necessary to know the lower limit of mass and force resolution above which SF and the stellar 
disk properties converge. Correctly determining the amount and location of SF is vital in 
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simulations in galaxies because SF 1) produces and distributes metals throughout the galaxy, 
2) affects the distribution of matter in the galaxy through feedback, and 3) enables us to relate 
simulations to observations. Therefore, we ask: what resolutions are sufficient to establish 
the galactic SF rates, the history of SF, the amount of stellar feedback, and the distribution 
of stellar and gaseous matter? To this end, we simulated a series of isolated galaxies of 
different masses at different mass and force resolutions and analyzed properties relating 
to SF. Isolated galaxies have the advantage of requiring less computational expense while 
allowing us to separate out environmental effects from those caused by different resolutions. 
We then analyze the SF and stellar feedback in the models and relate them to the history 
and structure of the galaxies to determine the resolution necessary for convergence and to 
describe the effects of resolution. 

The outline of our paper is as follows. In §2, we introduce the models and describe 
the computational methods used. The results of these simulations are described in terms 
of global SF (§3.1), stellar feedback (§3.2) and stellar distribution (§3.3). We address the 
physical connection between the SF, feedback, and structure and discuss implications of these 
results for large volume simulations in §4. We conclude §5 with a list of recommendations 
for future cosmological simulations. 



2. Methods and Models 

Our set of isolated galaxies ranges in mass from lO^M© to 10^'^ Mq. We simulated each 
galaxy at five mass resolutions with a range of 50 to 10^ DM particles and with the same 
number of initial gas particles. For simplicity, we refer to the mass resolution of a given 
simulation by the number of DM particles with the understanding that this also denotes 
the number of gas particles initially in the simulation. Based on our initial conditions, the 
mass of individual dark matter particles is 0.9 x Mh/Ndm whereas the mass of each gas 
particle is 0.1 x Mh/N^m, where Mh is the total halo mass and N^m is the number of dark 
matter particles. We chose to further focus on the IO^^Mq and the 10^°Mq galaxies. The 
former is representative of a Milky Way-mass galaxy. The latter is on the cusp of being able 
to form a disk through purely radiative cooling of primordial gas and, as will be shown, is 
particularly sensitive to SF and stellar feedback. In addition to the 50 to 10^ DM particle 
runs, we simulated these two galaxies using 10^ DM particles to check convergence at very 
high mass resolution. To examine the effect of force resolution, we re-simulated the lO^^M© 
and the IO^^Mq simulations with 10^ DM particles and a range of softening lengths between 
2.0 X 10^^ and 2.5 x 10~^Rmr, where Rmr is the viral radius, defined by the radius where the 
enclosed density is 200 times the critical density. A complete list of the simulations can be 



found in Table 1. 



2.1. Isolated Galaxy Model 



The i nitial conditions for o ur isolated galaxies are identical to those used i nlKaufmann et al 



(120071) and lStinson et al.l (120071 ) . They consist of an equilibrium, live DM halo (JKazantzidis et al 
20041 ). with a concentration of c = 8, where c = Rmr/Rs and Rs is the scale radius of the 
NFW profile. In order to establish equilibrium, the DM halo extends beyond Rmr at a 
slightly lower mass resolution, such that 72% of the DM mass and 90% of the DM particles 
are within the virial radius. The DM velocity distribution was set using the Eddington 
inversion method of Kazantzidis et al. 2004. Once the DM halo was created, 10% of the 
mass of the DM halo was converted to gas and distributed evenly over the same number 
of gas particles as DM particles. The gas was given the same NFW density profile as the 
DM and the temperature of each particle was chosen to establish hydrostatic equilibrium 
prior to cooling. The gas cloud was then spun with a uniform rotational velocity such that 
A = {jgas\E^/'^\){GM^/'^) = 0.039, where jgas is the average specific angular velocity of the 
gas and E and M are the total energy and mass of the halo. By creating our isolated 
galaxies from collapsing clouds of gas such as these, we avoid requiring the galaxies to be in 
equilibrium and enable the study of both a sudden surge of SF during the initial collapse 
and sustained SF later on. 



2.2. Numerics 



Our simulations were performed using GASOLINE, a parallel tree +SPH code (jWadslev et al 



20041 ). which is an extension of the A^-body gravity code PKDGRAV (IStadelll200ll ). A soft- 
ening length, e, is included in the calculation of gravitational force to decrease relaxation 
from close encounters of particles an d to make the simulation numerically feasible. We use 
spline kernel gravitational softening ( jHernquist fc Katzl 119891 ). The force on each particle 
is integrated using a timestep Atgrav = 'U\f^i where the timestep criterion rj = 0.175 and 
g is the acceleration . For gas particles, the timestep also satisfies the Courant condition 
( jWadsley et al.l 120041 1. Each gas particle has an SPH smoothing radius, h, over which gas 
properties are sampled. The value for h is set for each gas particle such that above a mini- 
mum value oi h cc e, the same number of particles are contained in each smoothing kernel; 
when h is equal to the minimum value, more particles may be contained within the kernel. 
In our simulations h > O.le and 32 particles are contained in each kernel when h > O.le. 
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Table 1. Simulations used in our analysis 

Mass [Mq] Softening Length [Rvir] Softening Length [pc] 

10^ 1.0 X 10-3 21 

10^°* 1.0 X 10-3 44 

10^° 2.5 X 10"^ 11 

10^° 2.5 X 10-3 110 

10^° 5.0 X 10-3 220 

10^° 1.0 X 10-2 440 

10^° 2.0 X 10-2 1760 

10^^ 1.0 X 10-3 96 

10^2* 1.0 X 10-3 206 

10^2 2.5 X 10-^ 51.5 

10^2 2.5 X 10-3 515 

10^2 5.0 X 10-3 1030 

10^2 1.0 X 10-2 2060 

10^2 2.0 X 10-2 8240 

10^3 1.0 X 10-3 444 



Note. — We simulated each of the listed galaxies with 50, 
100, 103, 10'^ and 10^ DM particles and an equal number of ini- 
tial gas particles. Galaxies marked with an asterisk were addi- 
tionally simulated with 10^ DM particles and an equal number 
of initial gas particles. For both the IQ^'^Mq and the lO^^M© 
galaxies, we additionally ran a suite of simulations with vari- 
ous softening lengths. The simulations with softening lengths 
not equal to 1.0 x 10-3/?^^^ are distinguished by the indented 
columns. 
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We implement gas cooling assuming ionization equilibrium and an ideal gas of primordial 
compositio n. The cooling ra te and the ion abundances are solved for using coUisional ioniza - 
tion rates (jAbel et al.l 119971 ). radi ative recom bination (JBlacklll98ll : IVerner &: Ferlandlll996l ). 
bremsstrahlung, and line cooling (jCenlll992l ). As in S06, the effects of the UV background 
radiation, metal line cooling and molecular hydrogen cooling are not included. 

We modeled SF using the sub-grid recipe outlined in S06. The formation of each star 
particle (representing a simple stellar population of a set metallicity) is a statistical process 
based on the properties of the local gas. In this recipe, gas particles must meet the following 
set of criteria to be capable of forming stars: the gas density must be greater than a minimum 
density, rimin, the temperature must be less than a maximurn temperature, T < Tmax = 
15, OOOK, and the gas must be in a converging flow (lKatall992l ). Whether a potentially star 
forming gas particle spawns a star particle in a time step. At, is determined stochastically 
with the following probability p: 



m, 



P 



gas 



^TT'star 



-C*At/tfor 



') 



(1) 



in which nigas is the mass of the gas particle, rristar is the initial mass of the potential star 
particle, c* is a constant star forming efficiency factor (the likelihood a valid gas particle 
will make a star particle in a formation time), and the formation time, tform, is equal to 
the dynamical time (S06). We chose rristar to be 40% of the initial gas particle mass as 
a compromise between the resolution of the galaxy's stellar component and computational 
expense. This set of minimum star forming crite ria and stochasti c SF p arameters leaves 
two tunable parameters: rimin and c*. The tests in lGovernato et al.l (|2007l) det ermined best 
values of rimtn = O.lcm"^ and c* = 0.05 and we use a iKroupa et al.l (119931 ) initial mass 
function to our star particles, as they did. 

Feedback from star particles represents the energy injected into the interstellar medium 
(ISM) by Type la and Type II supernovae. To calculate the ef fect of the released energy , 
we use the blastwave recipe described in S06 based on work of iMcKee fc Ostrikerl (119771 ). 
Energy from a given star particle is equal to the number of supernovae that would occur 
in a simple stellar population of the same mass and metallicity multiplied by a constant 
value, EsN = 4 x 10^° ergs per supernova (S06). The energy is then distributed through 
the surrounding gas particles and the cooling is disabled for a length of time and distance 
determined by the blastwave model. Feedback from stellar particles also returns mass and 
metals to the ISM. 
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3. Results 

3.1. Global Star Formation 

The global SFR of a galaxy (the mass of stellar particles formed in one year in the entire 
galaxy) is important to establish in part becaus e it can be asce r tained even in relatively 



low-resolution observational surveys of galaxies ( Bell et al.l l2003t iBrinchmann et al 



and from observations of high-redshift galaxies ( iMadau et al.l 1 19981 : iPapovich et al. 



20041) 



2001 



Shapley et al.ll200ll : IStanway et al.ll2003[ ). Furthermore, the global SFR is key to the integrity 
of the simulation because of the significant impact of stellar feedback on the structure of the 
galaxy. Through the injection of metals and energy by stellar feedback, gas content, chemical 
evolution, morphology, and future SFR of the galaxy are all dependent on the global SFR. 

Trends in the global SFR as a function of mass and force resolution can set lower 
limits on the appropriate resolution for simulations of different mass galaxies. The top panel 
of Figure [1] shows the average SFR over 3 Gyr for different masses of galaxies over our 
range of mass resolution. For galaxies with masses of lO^^M© or more, SFRs in simulations 
do not converge until 10^ DM particles are used; galaxies with these masses and 10^ DM 
particles or fewer do not form stars at all. For 10^^ Mq galaxies and smaller, low mass 
resolution simulations are still able to form stars but the SFRs differ greatly from those of 
simulations with high mass resolution. These low mass, low resolution galaxies form stars 
in our simulations because the absolute mass resolution is still relatively high, allowing the 
gas to more easily reach sufficient density for our SF criteria. For example, during the 
initial collapse about 20% of the gas particles in the 10^^ Mq, 100 DM particle simulation 
are sufficiently dense to form stars. In comparison, none of the gas in the lO^^M©, 100 DM 
particle simulation is capable of SF. In all cases, the SFR increases dramatically once 1000 
DM particles are used as more of the gas particles in all the simulations are able to reach the 
density threshold. The IO^^Mq galaxy is unique in showing a comparatively large increase 
in the SFR from 1000 to 10^ DM particles. This increase in SFR corresponds to an increase 
in disk strength and will be discussed in more detail in §3.3. As the number of initial gas 
particles is further increased to 10^, the SFR in this galaxy again declines. We will show 
that this can be linked to an increase in the effectiveness of stellar feedback (§3.2). 

The bottom panels of Figure [1] show the change in the average SFR over 3 Gyr for 
the IO^^Mq and lO^^M© galaxies with different numbers of particles and force resolutions. 
For the simulations with the highest mass resolution, values of e < 10~'^ R^ir result in very 
similar SFRs. For simulations with 1000 DM particles or fewer, however, the softening length 
begins to affect the SFR. For these simulations, softening lengths at both extremes result in 
smaller SFRs. Large values of e reduce the SF by spreading local density enhancements over 
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larger volumes. While at late times this is not sufficient to limit the amount of SF, during 
the initial collapse it prevents the formation of a high density region in the center of the 
galaxy. Because of this, about 20% fewer gas particles are able to reach SF densities in the 
e = 2 X 10~'^Rmr than in the e = 2.5 x 10~^Rmr simulation for both mass es of galaxies. Small 



value s of e, in contrast, are known to result in greater two-body heating ( ISteinmetz fc White 



19971 ). In our low mass resolution simulations (1000 DM particle or less) this has the effect 
of heating the gas and lowering the density such that fewer particles are capable of forming 
stars. As the softening lengths decrease from 10^^ to 2.5 x 10~^-Rj,jr, the fraction of gas that 
has sufficiently low temperature and high density to form stars after 3 Gyrs decreases from 
17% to 0% in the 10^^ Mq galaxy and from 11% to 2% in the lO^^M© galaxy. The global 
SFR of our simulated galaxies begins to converge for simulations with more than 10^ DM 
particles and for softening lengths smaller than 2.5 x 10~^Rmr- 

Star formation histories (SFHs) show the changes in the global SFR as the galaxy evolves 
and examining them gives us an idea of the global SFR on small time scale s. These may be ob 



servationally determined through stellar populations synthesis modeling (JBruzual fc Chariot 



20031 ) and are important for establishing the history of a galaxy and relating the SF to other 
process such as mergers, gas infall, disk instabilities, and stellar feedback. In our analysis, 
we calculate the global SFR for every 5 x 10^ years and normalize by the initial baryonic 
mass. Figure [2] shows both the effect of mass and force resolution on SFHs for the lO^^M© 
and the IO^^Mq galaxies. For the IO^^Mq galaxy, the SFH converge s for simula t ions w ith 



10^ DM particles or more. As seen in the cosmological simulations of iNaab et al.l (120071 ). in 
lower mass resolution simulations the initial burst of SF is delayed. SFHs of both galaxies 
converge at a softening length of approximately e = 10~^ Rmr- As was seen with the global 
SF, larger values of e reduce the SFR, particularly affecting the initial burst of SF because 
the sudden increase in density at the core of the collapsing cloud is effectively spread over a 
larger volume and shallower potential. 

The 10^^ Mq galaxy is extremely sensitive to changes in mass resolution. The SFH of the 
10^ DM particle simulation is dramatically different from all of the others and consists of a 
series of small bursts of SF. The 10^ and 10^ DM particle simulations both show a sustained 
burst of SF of approximately the same strength but offset in time by approximately 1.5 
Gyrs. Increasing the resolution from 10^ to 10^ DM particles decreases the strength of the 
burst as stellar feedback becomes more effective and the SFH qualitatively approaches the 
behavior of the 10^'^Mq galaxy. These sudden changes in the SFH leads to the differences 
in global SFRs seen in Figure [H The episodic SF of the 1000 DM particle simulation may 
be because this galaxy does not form a disk when simulated at low resolution, as will be 
discussed in §3.3. Instead, at low resolution this galaxy is supported in an elliptical shape 
by stellar feedback (§3.2), and the resulting SF shows a "breathing" pattern similar to that 
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seen in IStinson et al.l ( 120071 ) . 

In simulations with low mass resolution, there is a point where the large particle mass 

^0 



exceeds the Jeans masses of the gas particles capable of SF. In the case of the IO^^Mq 
galaxy, this happens for about 10^ particles or fewer, and in the IO^^Mq galaxy it happens 
for about 10^ particles and fewer. While the shift from bursty to sustained SF in the former 
case could be related to resolving the Jeans masses, it would have to be specific to this 
mass of galaxy as there is no obvious difference in the SFH for the latter case. To address 
the concern that as the mass resolution increased, the higher densities of the particles and 
lower temperatures would result in Jeans masses for our particles that were too small to be 
resolved, we further examined the Jeans mass of each particle. We found that in both the 
IO^^Mq, 10^ DM particle galaxy simulation and in the 10^ DM particle one, there was no 
difficulty in resolving the Jeans mass. For these two simulations, the lack of low temperature 
cooling kept the Jeans mass relatively high implying that the difference in SFHs could not 
be an effect of changing Jeans masses. 



3.2. Stellar Feedback 

The global SFR affects the structure of a galaxy through stellar feedback as energy 
from supernovae is injected into the surrounding media, heating the gas and, in some cases, 
expelling it from the galaxy. Not only must the SFR converge but the effects of SF through 
feedback must also conver ge. To determine the effect of feedback on gas in the galaxies, the 



gas was categorized as in iMac Low fc Ferraral (119991 ) according to its fate. Figure [3] shows 



the fractions of baryonic mass in the form of stars, accreting gas (cool, dense gas moving 
toward the center of the galaxy), and "blowaway" gas (unbound gas). The bulk of the 
remaining gas resides in the disk of the galaxy. In order to minimize any time dependency, 
the data was gathered after 2, 2.5 and 3 Gyrs and the fractions were averaged together. 
For both masses of galaxies, softening lengths of 5 x lO^^R^i^ are sufficient for convergence, 
representing the relative unimportance of spacial resolution on the efficiency of feedback 
above a certain threshold. Mass resolution, on the other hand, has a strong effect. As the 
number of particles increases, plumes of gas are more able to escape the disk. The fractions 
of gas converge for the 1000 DM particle simulation of the lO^^M© galaxy. 

For the 10^° M© galaxy, however, the mean SFR does not converge even with 10^ DM 
particles. For this galaxy, the effect of feedback increases, causing much more gas to be 
expelled from the galaxy and less to be accreted. Figure H] shows the large plume of gas 
driven from the 10^ DM particle galaxy. Compared to the lO^^M© galaxy simulated with 
10® DM particles, this plume is a much larger fraction of the total mass of the galaxy. 
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As is expected from the comparatively shallow potential, gas is much more easily expelled 
from this galaxy. This shallow potential, in combination with the changing morphology of 
the disk (§3.3), make feedback in the lO^'^ Mq galaxy particularly sensitive to resolution. 
As the number of DM particles increases from 10^ to 10^, the mass of gas expelled sharply 
increases. However, the 10^ DM particle simulation has less total SF than the 10^ DM particle 
simulation at almost all times (Figure [2]), indicating that the large amounts of feedback in 
the former are not caused by greater total amounts of SF. Indeed, feedback appears to be 
suppressing SF in this galaxy. Instead, we suggest that the different amounts of expelled gas 
are because of changes to the gas disk. Increased mass resolution results in thinner disks, as 
will be shown in the following section. Additionally, using more particles results in smaller 
average smoothing lengths causing decreased amounts of numerical viscosity in high mass 
resolution simulations. The exact effects of lowering the numerical viscosity in simulations 
are poorly understood, making it difficult to disentangle these two potential causes. Both 
thinner disks and lower numerical viscosity, however, would be expected to aid the expulsion 
gas by reducing the force necessary for the winds to push out of the disk. 



3.3. Disk Structure 

While the global SF may be sufficient information to analyze large-scale simulations 
of structure formation, when studying individual galaxies the spatial distribution of SF is 
important. The distribution of stars and overall morphology of galaxies is ob servationally 



used to characterize galactic type, indicates the merger history of the galaxy (ICouch et al. 



19981 ). and is used to establish trends such as the Tully-Fisher relationship. We examine the 
effects of mass and force resolution on the location of SF by looking at the properties of the 
stellar disks. In particular, we relate the mass and force resolution to the resulting rotation 
curves, scale heights, and ellipticities of our galaxies. 

Figure 5 shows the effect of resolution on the rotation curve. The circular velocity, 
Vc, is defined such that Vc = \/GM/r. Almost of our galaxies have a high central mass 
concentration as represented by the steep increase in Vc at small radii. The c oncentrations are 



unaffe cted by increasing mass resolution. This contrasts with the findings of iGovernato et al. 



teOOTl ). in which an increase in resolution similar to the increase from 10^ to 10^ particles in 
our IO^'^Mq simulation resulted in a significant decrease in the inner circular velocity. The 
difference likely results from the isolated environment of our simulations during which gas 
collapses into pre-existing DM halos and within which there are no outside influences such 
as mergers or tidal torques to affect the galaxy's angular momentum. 

One simulation that does show a signiflcant decrease in the inner circular velocity is 
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the 10^° Mq galaxy simulated with DM 10^ particles. In this galaxy, the simulation with 
10® DM particles is strikingly less centrally concentrated than those of lower resolution. 
The low concentration is the result of the highly effective stellar feedback in this particular 
galaxy, as seen in t he la rge p lume of escaping gas iii Figure HI Similarly to the galaxies in 
Mashchenko et al.l (120081 ) and iGovernato et al.l (J2010l ). the large amounts of stellar feedback 
appear to have reduced and flattened the central dark matter concentration. Within the 
central 30 to 300 pc of the lO^^M© galaxy, the logarithmic slope of the dark halo decreases 
from —1.01 in the 10^ DM particle simulation to —0.51 in the 10® DM particle one. A 
second key effect of the strong supernova-genera ted outflows, whi c h is al so seen in the high SF 
threshold cosmological simulations presented in IGovernato et al.l (120101 ) , is the extremely low 
SF efficiency at the center of the galaxies that pre vents the formation of a s tellar cusp. Both 
the cosm ological simulations of d warf galaxies by iMashchenko et al.l (120081 ) (2 x 10^ Mq at z 
= 5) and IGovernato et al.l (120101 ) (3.5 x IO^^Mq at z = 0) included gas cooling through metal 
lines and a high density threshold for star formation {rjmin = 10 amu cm~^ and rjmin = 100 
amu cm~^, respectively), in addition to exceptionally high mass resolution. Despite the fact 
that our simulations allowed SF in lower density gas and not did not include a prescription 
for cooling below 10'^ i^, this galaxy demonstrates how feedback can affect the DM as well 
as the baryonic distribution of matter in isolated simulations. 

In general, force resolution affects the rotation curves more than mass resolution. While 
low force resolution simulations result in rotation curves more similar to observed galaxies, 
this is a numerical artifact. The large softening lengths simply cause the gas properties to 
be averaged over larger volumes and prevent a concentrated central density from forming. 

As shown in Figure [6l both mass and force resolution affect the scale height of the stellar 
disk. A decrease in the mass resolution causes the scale heights of the 10^'^Mq and IO^^Mq 
galaxies to be larger at most radii, indicating puffier disks. As the number of particles 
increases, two-body heating decreases, and thinner disks are able to form. Varying the 
softening length has a negligible effect on the scale height of the disk of the lO^^M© galaxy. 
For the lO^^M© galaxy, on the other hand, lower force resolution has a similar effect to 
lower mass resolution: it increases the scale height. The puffier disks at low force resolution 
are produced as the shallower potential allows the feedback to become more effective. This 
appears to be another situation in which the IO^^Mq galaxy is very sensitive to numerical 
parameters and is connected to the wide range of feedback efficiencies seen in this mass of 
galaxy. 

In order to combine the radial stellar disk profile and the scale height of the disk into 
a single parameter representing disk strength, we calculate the ellipticity of the stellar dis- 
tribution as viewed edge on. To do this, we chose the major axes to be r^go, the radius 
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at which 90% of the stellar mass is enclosed, and the minor axes to be 2;*,9o, the height at 
which 90% of the stellar mass is enclosed. Using r* go and 2*90, we defined the ellipticity 



to be e 



«,90 



where e = 1 is a disk and e = is a sphere. Using the ellipticity 



of the galaxy, measured in this way, as a proxy for the strength of the disk is an ad hoc 
method and was chosen simply as a straight forward way of quantifying the galaxy's appear- 
ance. The upper panel of Figure [7] shows the ellipticity of galaxies of different masses and 
mass resolutions. Across all resolutions, the IO^Mq galaxy has a more spherical shape. The 
IO^^Mq galaxy is the lowest mass galaxy to form a disk, the formation of which is strongly 
dependent on the number of particles in the simulation. In general, the strength of the disk 
(i.e. greater ellipticity) increases sharply with mass resolution. For the lO^^M© g^l^^Y) the 
ellipticity also roughly increases with force resolution (Figure [7]) as the more concentrated 
potential causes the disk scale height to decrease. The ellipticities of the IO^^Mq galaxy 
simulated with 10^ and 10^ DM particles, however, turn over at the highest force resolution 
(Figure [7]). This is the result of the decreasing disk scale lengths at high force resolutions 
while the scale heights remains largely unaffected by changes to the force resolution. The 
disk parameters, including rotational curves, scale height, and ellipticity, converge for the 
10^^ Mq and 10^'^ Mq galaxies for softening lengths between 10~^ and 10~^ Rmr- Generally, 
disk parameters converge in simulations of 10^ DM particles and higher. 



4. Discussion 



In our simulations, we demonstrate a close relationship between resolution, stellar feed- 
back, spatial distribution of star formation, and the am ount of star formation. Cooling 
in combination with angular momentum produces disks (jOstriker fc Binneyl 119891 ). Angu- 
lar momentum is artificially reduced in simulations with few particles, shrinking th e size 
of disks formed if not preventing disks from forming entirely ( IGovernato et al.l 120041 ). The 
higher resolution models create the thinnest (i.e. highest-ellipticity) disks, and form the 
most stars due to the presence of the highest density gas. Increasing the resolution results 
in sustained, constant SF with a feedback-driven plume of gas. In high-mass galaxies, this 
occurs in simulations with 10'^ particles where the SF and the effects of feedback begin to 
converge. In low- mass, lower- resolution galaxies with thick disks (i.e. low-ellipticity), SF is 
episodic as feedback shuts off SF for short periods of time. For these galaxies, it is most likely 
a combination of the relatively large disk scale height, the more spherical structure, and the 
greater artificial viscosity that prevent the heated gas from puncturing through the galaxy 



in the form of a galactic wind ( iMac Low &: McCrayl Il988l : iMac Low et al 



stellar feedback causes a temporary halt to SF in the galaxy (IStinson et al. 



2003). 



1989). Instead, 
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The IO^^Mq galaxy is at the transition between these two behaviors and its shallow 
potential makes it particularly sensitive to stellar feedback. Increasing the mass resolution 
from 1000 to 10^ DM particles results in the formation of a disk. Further increasing the 
number of DM and initial gas particles from 10^ to 10^ produces increased gas loss through 
stellar feedback. At our highest mass resolution, feedback is efficient enough to decrease the 
central concentration of the halo and to produce a slowly rising rotation curve. Over the 
same range of mass resolution, the SFH goes from small episodic bursts, to one extremely 
large burst, to a smaller burst followed by sustained SF. These changes appear to be linked; 
low resolution causes the galaxy to be less disky, making supernova feedback expand the gas 
spherically outwards and ending SF until the gas is able to one again cool and compact. At 
the highest resolution, the galaxy contains a strong disk and outflows of gas concentrated 
along the poles. While these outflows do not cause a galaxy-wide halt to SF, they do lower 
the concentration at the center of the galaxy. We expect that it is the combination of smaller 
disk scale-height and lower artiflcial viscosity in the higher mass resolution simulations that 
promote the formation of the outflows. The force resolution also has relatively strong effects 
on this galaxy, as large values of e cause thicker disks. In fact, even at our highest resolu- 
tion, we were unable to establish the convergence of the SFR, SFH, or the fate of the gas. 
This sen sitivity extends eve n beyond resolu tion: changing the i nitial mass function from 
Kroupa (JKroupa et al.lll993l ) to Miller-Scalo (JMiller fc Scaldll979l ) produced similar changes 
to the SFH. The susceptibility of simulations of this mass of galaxy to small changes in the 
parameters is at least in part due to its low virial temperature; its virial temperature of 
approximately 2.5 x 10^ K is close to lO^K, which is the temperature regime at which metal 
line cooling becomes important. The set of simulations examined in this paper does not 
include cooling from metals or H2 and low temperature cooling is necessary to study these 
low-mass galaxies in greater depth. 

There is some evidence that the strong sensitivity of this r riass of galaxy i s not unique 
to our simulations but may be based on a physical reality. IStrigari et al.l ( l2008l ) found 
that the Milky Way satellites share a common mass of lO^M© within their central 300 
parsecs despite covering a IO^Lq range in luminosity. Simulations predict that each of these 
satellites originated from galaxies with a tota l mass prior to interacting with the Milky Way 
of approximately lO^M© (IStrigari et al.ll2008l ). Despite the similar halo mass, differences in 
history and environment were capable of causing extremely different SFHs. Although we are 
considering galaxies a factor of 10 more massive, it is difficult to match the masses, both 
because of the extrapolation of the total halo mass from the observed mass within 300 pc and 
the extrapolation from our isolated simulations to a satellite halo. Therefore, the observed 
IO^Mq galaxies could be in the same precarious position as our IO^^Mq galaxies. If such 
a wide range of luminosities in the Milky Way satellites can be produced by similar initial 
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mass galaxies undergoing different histories, it is possible that the sensitivity of our 10^^ Mq 
reflects the sensitivity of real galaxies of similar mass. Our studies demonstrates that slightly 
different amounts of feedback and ease of mass loss can completely change the structure and 
SFH of this mass of galaxy. While in this study, those differences were the result of changes 
to resolution, these could also be caused by differing merger and gas accretion histories. 

One question when examining these results is how universal the requirements for con- 
vergence are. In particular, will the same resolution be necessary in galaxies created from 
different initial conditions, galaxies using different star formation parameters, and galaxies 
that are evolved with different feedback recipes? Similar mass resolution was required for 
convergence of the total amount of SF in cosmological runs of Milky Way mass galaxies 
(JBrooks et al.l 120071 : iGovernato et al.ll2007l ). with slightly different feedback energies and SF 
efficiencies, implying that the resolution necessary for convergence is mostly independent 
of the initial conditions and of small changes in those two SF parameters. Feedback reg- 
ulated SF tends to reach an equilibrium with increased SF efficiency resulting in stronger 
feedback which then hinders the formation of more stars. Based on the relative insensitivity 
to these parameters, we can assume that similar mass resolution would be necessary even 
with different values for these parameters. The minimum density for SF, rjfnin, is the single 
parameter one would expect to be strongly linked to resolution. As lower mass and spatial 
resolution results in less high density gas, increasing rjmin requires resolving the structure of 
the ISM to a greater degree (IGovernato et al.ll2010l ). Simulations with high values of rjmin, 
therefore, should need greater resolution for convergence. A more cor nplete discus s ion of the 
relationship between efficiency, rjmm and resolution may be found in ISaitoh et al.l ( 120081 ). 



A more difficult issue is the robustness of our results across different feedback recipes. 
Although most low-resolution SPH simulations tend to have depressed amounts of SF caused 
by the presence of a smooth, thick gas disk, the actual mass resolution necessary for con- 
vergence appears to be high ly dependent on the feedback recipe. For example, in the 
Springel fc HernquistI (l2003al ) recipe used in GADGET individual gas particles represent a 
multi-phase ISM and some of the energy from feedback is transferred to selected nearby gas 



particles in the form of kicks to their velocities. As seen in Figure 12 of ISpringel &: Hernquist 



( l2003al ). this approach has the advantage of being relatively insensitive to mass resolution. 
By imposing an equation of state upon the gas particles, self-regulation of SF is assumed; 
by allowing sub-grid multi-phases of the ISM within gas particles, different phases of the 
ISM need not be resolved. The cost of these advantages, however, is ease in describing non- 
equihbrium feedback configurations. For SF to be shut off or for superbubbles to be formed, 
a model for galactic winds must be assumed. One approach when mod eling galactic winds is 



to tera porarily decouple the winds hydrodynamically from the galaxy (ISpringel &: Hernquist 



2003aJ). While this makes the simulations largely insensitive to resolution, the results are 
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very sensitive to the model parameters. In models in which the wind does remain coupled to 
the gas, decreasing resolution induces greater mass loss in the win ds and, in more massive 
galax ies, causes the winds to have less effect on the disk structure (JDalla Vecchia fc Schaye 
20081 ). Despite the differences in feedback recipes, our results agree qualitatively in that 
higher mass and spatial resolution results in stronger d isk structure and (aside frorn the 
10^*^ M0, 10^ DM particle halo) higher SF. However, while ISpringel fc HernquistI (J2003al ) was 
able to establish convergence of the SFR of the IO^^Mq halo, our simulations did not. We 
believe this is because our simulations of this galaxy are in highly non-equilibrium states 
between star formation and supernova feedback, whereas some degree of self-regulated star 
formation is assumed in the ISpringel fc HernquistI ( l2003al ) simulations. We find in general 
that for cases other than the 10^^ Mq halo, we need at least a factor four higher mass res- 



olution in our simulations for convergence of the global SFR than in ISpringel fc Hernquist 



The effect of resolution on global SF and disk structure has implications when studying 
SF in cosmological simulations. Cosmological simulations generally use the same mass of 
particles and the same spatial resolution throughout. This causes less massive halos to be 
made-up of relatively few particles, each of whose softening lengths are a larger fraction of the 
virial radi us, and results in less accurate SF for these halos. Several cosmological simulations 
using the ISpringel fc HernquistI (l2003af) feedback recipe or modifications thereof have been 



used to study the cosmological SFH (ISpringel &: HernquistI l2003bl : IChoi fc Nagamind 12009 



Schaye et al.ll2009l ). In these studies, the authors balance the need to simulate a large volume 
at z = with the importance of high mass resolution at high redshift by stitching together a 
series of simulations of increasing size and decreasing mass resolution at lower redshifts. Over 
each redshift interval, the cos mological SFRs within simula ti ons of different mass resolutions 
are checked for convergence. ISpringel fc HernquistI ( l2003bl ): ISchaye et al.l (120091 ) noted that 
while low resolution results in a decrease in the SFR at high redshifts, the amount of free 
gas left over may result in artificially high SFRs at z=0. 

One example of a large , cosmological SPH simulation with similar SF an d feedback 



Keres et al.l ( 120051 ). In their primary simulation, iKeres et al.l ( 120051 ) simulated 



recipes is 

galaxy formation with in a 22.222 



1 ^ M pc comoving cube using the star formation and 



feedback recipes from iKatz et al.l (Il996[ ). Stars formed probabilistically from eligible gas 
particles and supernova feedback energy is transferred to the surrounding gas particles in the 
form of heat. This feedback recipe is similar to ours in that it took place over single ISM phase 
gas particles and did not assume regulated SF, but the lack of the blastwave formulation 
would have c aused the fe e dback -heated gas particles to cool more quickly. According to our 
tests, in the iKeres et al.l ( 120051 ) primary cosmological simulation (L22/128), galaxies with 
masses of IO^^Mq and above were simulated with a sufficient number of particles to judge 
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the global SFR. For less massive galaxies, however, the mass resolution was low enough to 



throw into question their SFR. In their simulation, galaxies with masses close to IO^^Mq 
were simulated with roughly 100 particles. Our isolated galaxies of similar mass and particle 
number have 100 times less SF than their high-resolution counterparts over three billion 
years, although the difference decreases after the initial collapse and burst of SF. We estimate 
the effect of using low mass resolution on low mass galaxies when calculating the total SFR 
in a simulation by multiplying the equilibrium SFR in each indivi dual mass galaxy by the 
number of those galaxies expected from the galactic mass function (JReed et al.ll2007l ). When 
comparing the total SFR after the initial collapse (from 1.5 to 3 billion years) generated 
by our best resolved galaxies to that generated by galaxies with particle masses of IO^Mq 
(close to L22/128 resolution), we find our well-resolved simulations to produce about 10% 
more SF. The net difference in SF is significantly smaller than the discrepancy for low mass 
galaxies because most stars form in L^ galaxies, which were simulated at higher resolution. 
It should be noted, however, that at higher redshift, low- mass galaxies make up a higher 
percentage of the overall mass and the fraction of SF that is poorly constrained would 
increase. Furthermore, the relationship between resolution and the global SFR, particularly 
in the low-mass galaxies, is not a linear one and intermediate resolution simulations could 
actually produce greater discrepancies in total SFRs. 

Because of the difficulties in comparing simulations of isolated galaxies to cosmological 
simulations as well as the difference between the feedback recipes used in the simulations, 
this study d o es no t ser ve as a direct esti r nate o f the possible error in the SFR. As seen in 
Mayer et al.l ( 120081 ) and iGovernato et al.l (120081 ) . excess angular momentum loss may occur 
at different resolutions in cosmological simulations versus isolated galaxies. Without the ef- 
fects of large-scale tidal torques from cosmological simulations, our isolated galaxies may not 
have the same angular moment um, morphology, o r disk structure. Furthermore, the faster 
radiation of feedback energy in iKeres et al.l (120051 ) would result in decreased regulation of 
SF, particularly at the low mass end, and their feedback recipe may be affected differently 
by resolution. This exercise does highlight, however, the importance of resolution in cosmo- 
logical simulations and illustrates how knowledge of the effects of resolution over a range of 
masses can be used to generate approximate limits for the analysis of large scale simulations. 



Conclusions 



With limited computational resources, a compromise is always necessary between sam- 
pling a large number of galaxies and resolving those galaxies well. The practical question 
is not at what level of resolution do all measurements of the galaxy converge, but what 
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minimuin resolution is needed to answer a specific question. Studying the bright end of the 
galactic luminosity function, for example, requires very different resolution than studying 
the SFHs of dwarf galaxies. 

Based on our isolated galaxies, we recommend resolutions for various modeling goals in 
cosmological SPH simulations for the S06 SF and feedback recipes. These recommendations 
may also be applicable to any SPH simulation with single-phase gas particles in which 
supernova feedback energy is transferred to the surrounding gas particles. As these are based 
on isolated galaxies, they come with a number of caveats when extrapolating to cosmological 
simulations. Namely, our galaxies have purely spherical accretion, have no mergers, and have 
no outside tidal forces. These limitations affect the properties of accreted gas as well as the 
amount of angular momentum transferred to and from the galaxy. Our initial conditions 
also only sample one standard angula r momentum (A = 0.039) and one mass profile (NFW) 



and assume the iKroupa et al.l (119931 ) initial mass function. Despite these caveats, our set 
of simulations covers a wide range of masses, mass resolutions, and force resolution and 
produce both "bursting" and sustained SF at different stages. Together they provide a point 
of reference for determining the necessary resolution of galaxies. With these limitations in 
mind, we make the following recommendations. 

• For simulations similar to ours with gas particle masses less than IO'^Mq, the Jeans 
mass of star forming gas is resolved This implies that our values of Umin and c* in the 
SF recipe are suitable for our range of resolutions using our cooling recipe. 

• Fewer than 10^ DM particles or initial gas particles or softening and smoothing lengths 
greater than 10~^Rmr do not produce accurate global SFRs and supernova feedback. 
Over a range of masses, we find that the average global SFR converges for isolated 
galaxy simulations of 10^ particles and for softening lengths smaller than 10~^Rmr- 
Gas loss and accretion also converge for similar values (except in the 10^ DM particle, 
10^^ Mq simulation). While the connection between mass loss and SF is clear, it is 
remarkable that the mass loss convergence happens at the same values over almost the 
entire range of galactic masses studied, despite the differences in the galactic potentials. 

• SFHs for galaxies do not converge for simulation using less than 10^ DM particles 
and less than 10^ initial gas particles or for simulations with softening lengths larger 
than 10~^ Ryir- In galaxies near IO^'^Mq, the interaction between feedback and the 
shallow potential wells create spurious effects and prevents convergence at all mass 
resolutions. In general, increasing mass resolution results in more peaked bursts of 
SF followed by a smooth decrease in the SFR. The IO^^Mq galaxy, however, shows 
a wide range of behavior over the mass resolutions tested with no clear convergence. 
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For this mass of galaxy, raising the number of particles causes the creation of a disk 
and raising it further results in a sudden increase in the amount of mass lost through 
supernova feedback. The SFHs of this particularly sensitive galaxy reflect this: at low 
mass resolution the lack of a coherent disk results in episodic SF; slightly better mass 
resolution simulations show higher peaked SFRs; the highest resolution galaxy shows 
less SF as stellar feedback becomes more efficient. 

The spatial distribution of stars is not robust for softening lengths of e > 10~'^Rmr or 
in simulations with fewer than 10^ particles. Low mass and force resolution affect the 
spatial distribution of the stars by making the galaxy less disky with a less steep central 
density profile. Improving the mass resolution reduces the angular momentum loss and 
shrinking the softening length similarly reduces artificial pressure and allows disks to 
collapse more fully. Disk properties in galaxies with 10^ particles are only starting to 
converge and 10^ particles are generally necessary for the results to be robust. 



For galaxies much smaller than 10 Mq, low temperature cooling from metals (iMashchenko et al. 



20081 ) and H2 is necessary for gas to reach temperatures below the virial temperature. This 
is particularly important when predicting the effect of stellar feedback in dwarf galaxies or 
at high redshift. High resolution simulations in combination with low temperature cooling 
and H2 formation would better allow us to probe the star formation within dwarf galaxies, 
including the effects of feedback on rotation curves and episodic versus sustained star for- 
mation. In future work, we will examine the effect of H2 on cooling in galaxies of different 
masses and present an H2-based SF model 
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Fig. 1. — Average global SFR as a function of particle number and force resolution over 
3 Gyrs. The top panel shows the average global SFR for galaxies over a range of mass 
resolutions, normalized by the average global SFR of the 10'^ DM particle simulations for 
each mass. This normalization was chosen to highlight the effect of the changing mass 
resolution, particularly in lO^^M© galaxy. Each curve represents a different mass galaxy 
and the x-axis delineates the number of DM particles used. For these galaxies, the force 
resolution is kept at lQ~^Rmr- The lower two panels show the effects of both force resolution 
and the number of DM particles for the IO^^Mq (left) and lO^^M© (right) galaxies. In these 
panels, each line represents a different number of particles. 
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Fig. 2. — SFHs for lO^^M© and IO^^Mq galaxies for different force and mass resolutions. The 
SFHs are calculated using bins of 5 x 10^ years and are normalized by the initial baryonic 
mass (Mb). The two left panels show the effect of the number of DM particles on the SFR 
for a 10^'^ Mq (top) and 10^*^ Mq (bottom) galaxy with force resolution of 10'^^ Rmr- The two 
right-hand panels the show the effect of varying force resolution for a 10^'^Mq (top) and 
IO^^Mq (bottom) galaxy with 10^ DM particles. 
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Fig. 3. — Fate of the gas in 10^^ and lO^^M© galaxies for different force and mass number 
resolutions, with data collected at 2, 2.5 and 3 Gyrs and averaged together. The left panels 
are for a IO^^Mq (top) and lO^^M© (bottom) galaxy with force resolution of 10~^Rmr and 
show the effect of different numbers of particles. The two right-hand panels are also for a 
lO^^M© (top) and lO^^M© (bottom) galaxy but with 10^ DM particles and varying forc e 
resolution. The categorizations for the gas are adopted from iMac Low &: Ferraral (119991 ). 
The accreting gas is designated as cold, dense gas moving toward the center and the stellar 
fraction is the relative amount of gas that formed into stars. "Blowaway" denotes the fraction 
of baryonic mass in the form of gas permanently lost to the galaxy. 
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Fig. 4.— Distribution of gas in the IO^^Mq galaxy with 10^ DM particles (left), the 1O^°M0 
galaxy with 10^ DM particles (middle), and the 12^^ Mq galaxy with 10^ DM particles (right) 
after 3 Gyrs. The grayscale corresponds to surface density of gas for pixel on a logarithmic 
scale. Because of different virial radii, the same color in the 10^'^Mq galaxy represents a 
factor of 22 times more in total mass per pixel. Both simulations of the 10^°Mq galaxy show 
escaping plumes of gas caused by feedback but the 10^ DM particle simulation undergoes 
much greater mass loss. 
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Fig. 5. — Rotation curves of the 10^^ and IO^^Mq galaxies for different force and mass 
resolutions after 3 Gyrs. The two left panels show the effect of the number of particles on 
the profiles for a lO^^M© (top) and lO^^M© (bottom) galaxy with e = 10~^Rmr- The two 
right-hand panels show the same lO^^M© (top) and lO^^M© (bottom) galaxy but simulated 
with 10^ DM particles and a range of force resolutions. 
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Fig. 6. — Effect of different mass and force resolutions on the scale height of the stellar disk 
as a function of radius for the 10^^ and lO^^M© galaxies. The stellar scale height is defined 
as the height at which 90% of the stellar matter is enclosed and is plotted as a function of the 
radius for different resolutions and masses. The two left panels show the scale heights after 3 
Gyrs for a 1()^'^Mq (top) and 10^^ Mq (bottom) galaxy with force resolution of e = lO^^Rmr 
simulated with different numbers of DM particles. The two right-hand panels show the effect 
of varying the force resolution for a IO^^Mq (top) and lO^^M© (bottom) galaxy simulated 
with 10^ DM particles. 
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Fig. 7. — Ellipticity (e, defined for tlie stellar matter as in §3.3) of the galaxies viewed edge 
on after 3 Gyrs. In the upper panel, the ellipticity is plotted against the number of DM 
particles in the simulation for each mass of galaxy. The lower panels show the effect on the 
ellipticity of changing the spacial resolution for a lO^^M© (left) and a IO^^Mq (right) galaxy 
with 10^ DM particles. 
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